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With ultracold 87 Sr confined in a magic wavelength optical lattice, we present the most precise 
study (2.8 Hz statistical uncertainty) to-date of the 1 So - 3 Po optical clock transition with a detailed 
analysis of systematic shifts (20 Hz uncertainty) in the absolute frequency measurement of 429 228 
004 229 867 Hz. The high resolution permits an investigation of the optical lattice motional sideband 
structure. The local oscillator for this optical atomic clock is a stable diode laser with its Hz-level 
linewidth characterized across the optical spectrum using a femtosecond frequency comb. 



In recent years, optical atomic clocks have become 
increasingly competitive in performance with their mi- 
crowave counterparts. The overall accuracy of single 
trapped ion based optical standards closely approaches 
that of the state-of-the-art cesium fountain standards 
p], Q . Large ensembles of ultracold alkaline earth atoms 
have provided impressive clock stability for short aver- 
aging times, surpassing that of single-ion based systems. 
So far, interrogation of neutral atom based optical stan- 
dards has been carried out primarily in free space, un- 
avoidably including atomic motional effects that typically 
limit the overall system accuracy U S • An alternative 
approach is to explore the ultranarrow optical transitions 
of atoms held in an optical lattice 0, 0, Q • The atoms 
are tightly localized so that Doppler and photon-recoil 
related effects on the transition frequency are eliminated 
Meanwhile, the trapping potential is created by laser 
light at a carefully chosen wavelength (A mag i C ) such that 
it has essentially no effect on the internal clock transition 
frequency. Additionally, the increased atom-probe laser 
interaction time enabled by the lattice confinement per- 
mits a full utilization of the narrow natural linewidth. 
This optical lattice approach using neutral atoms may 
provide the best possible combination of clock stability 
and accuracy. Such a proposal has been under intensive 
investigation in the case of the doubly forbidden 1 So - 
3 Po transition in the fermionic Sr isotope, 87 Sr jtj IToj. 
Similar work in Yb is also in progress 11, 12 . 

The first high precision absolute frequency measure- 
ment was recently reported for the 87 Sr 1 6'o - 3 Po clock 
transition using a GPS-based frequency reference Q. 
However, establishing a new standard demands that dif- 
ferent groups study potential systematic errors. In this 
Letter, we present the most precise study to-date of this 
ultranarrow clock transition in a magic wavelength opti- 
cal lattice, with a direct reference to the NIST Fl Cs 
fountain clock. We have investigated systematic fre- 
quency shifts including those originating from atomic 
density, wavelength and intensity of the optical lattice, 
residual magnetic field, and probing laser intensity. We 
determined the absolute frequency of the 87 Sr clock tran- 
sition to be 429, 228, 004, 229, 867 ± 20 (sys) ± 2.8 (stat) 
Hz. Furthermore, we apply this optical clock resolution 
towards a study of ultracold atomic motion in an anhar- 



monic optical potential. 

To exploit the long coherence time provided by the ul- 
tranarrow (transition linewidth 7 ~ 1 mHz) 1 So - 3 Po 
clock transition, we have developed a cavity stabilized 
diode laser at 698 nm. The laser is locked to an iso- 
lated, passive reference cavity with finesse of ^250,000 
and cavity linewidth of ~20 kHz. This reference cav- 
ity is mounted vertically and supported symmetrically 
with respect to the cavity length, which minimizes the 
length sensitivity to environmental vibrations 13]. The 
entire system is relatively simple and compact, occupying 
less than 1 m 3 . To characterize this laser, we compare 
it to a highly stabilized NdiYAG CW laser (in an ad- 
jacent laboratory) which has consistently demonstrated 
a sub-Hz linewidth. This comparison was made possi- 
ble with a phase-stabilized femtosecond frequency comb. 
The repetition rate (100 MHz) of the octave spanning, 
self-referenced Ti: Sapphire fs comb was tightly locked 
to the cavity stabilized 698 nm diode laser, coherently 
transferring the diode laser stability to each of the ~10 6 
modes of the fs frequency comb. A heterodyne beat sig- 
nal between the sub-Hz Nd:YAG reference laser and a 
corresponding fs comb mode at 1064 nm, shown in Fig. 
1, revealed a linewidth of less than 5 Hz for the 698 nm 
diode laser and demonstrates optical coherence transfer 
at the 10 -14 level between remotely located lasers of dif- 
ferent colors. Additionally, the final stability of the 698 
nm laser system was confirmed to be limited by thermal- 
mechanical noise of the mirror substrates in the Zerodur 
reference cavity [l^ and will be improved with new mir- 
rors. The frequency drift caused by the material creep of 
the cavity spacer was carefully compensated via measure- 
ment of the optical frequency relative to a Cs fountain- 
calibrated hydrogen maser using the same fs comb. 

Preceding the 1 S'o - 3 Po spectroscopy, 87 Sr atoms un- 
dergo two stages of laser cooling (using the 1 So - 1 P\ and 
then the 1 So - 3 Pi transitions - see Fig. 2) to prepar e ap- 
proximately 10 6 atoms at fiK temperatures [lM Er l fl7| . 
Loading of atoms from the MOT to a 1-D optical lat- 
tice operating at A mag i C « 813 nm allows simultaneous 
cooling and trapping on the 1 5'o - 3 Pi transition so that 
spatial mode matching between the MOT and lattice is 
not critical. 10% of the atoms (> 10 5 ) are loaded from 
the red MOT to the lattice, which has a lifetime of 1 s. 
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FIG. 1: (color online) Measurement of probe laser (698 nr 
linewidth by comparison to a reference laser at 1064 nm ush 
an octave spanning fs laser. Linewidth of the heterodyne be 
at 1064 nm is 3 Hz. 




FIG. 2: (color online) Partial diagram of 87 Sr energy levels 
and the schematic of the optical layout delivering the probe 
and lattice laser to the atoms. P: polarizer, DM: dichroic 
mirror, PMT: photo-multiplier tube. The strong longitudinal 
and weak transverse confinement of the 1-D lattice creates 
2-D trapping regions. 



Inside the lattice, the axial trap oscillation frequency is 
Q = 2tt x 80 kH z, corresp onding to a Lamb-Dicke param- 
eter of r] = k p ^/h/ (2mVt) = 0.23. Here k p — 2tt/\ p , X p is 
the probe wavelength, 2ttH is Planck's constant, and m is 
the atomic mass. In this regime, spectroscopy of the clock 
transition can be performed nearly free of any Doppler or 
photon recoil shifts. To avoid residual transverse Doppler 
effects (transverse oscillation frequency is 500 Hz), the 
probing laser is carefully aligned to co-propagate with 
the lattice laser (Fig. 2), with both lasers being focused 
to a ~70 jam beam diameter at the trap. The probing 
laser is operated with optical powers near 5 nW. Spec- 
troscopy typically consists of a probing time of 10-40 ms 
followed by a 2 ms illumination of the atoms by 679 nm 
light, which pumps the excited state 3 Po population to 
the 3 Si state for shelving to the 3 P2 state after sponta- 
neous decay (Fig. 2). Because of spatial inhomogeneity 
in the Rabi excitation frequency, we repeat the iteration 
of probe and pump pulses approximately 20 times to en- 
hance the signal size. Finally, we detect the ground state 
1 6'o population by resonantly exciting the strong 1 S'o - 
1 Pi transition at 461 nm and counting fluorescence. 

The resolution of lattice spectroscopy enables us to in- 
vestigate the motional states of atomic Sr along the lon- 
gitudinal axis of the 1-D lattice. A deep lattice potential 
(where the populated motional quantum states are well 
below the energy barrier between adjacent sites) can be 
well approximated by harmonic confinement leading to 
equally spaced motional states. These motional states al- 
low inelastic excitations which are red and blue detuned 
from the elastic, purely electronic 1 5'o - 3 Po excitation, 
with amplitudes suppressed by ?7 2 /4 in the Lamb-Dicke 
regime. A measure of the amplitudes of the red and blue 
sidebands directly gives the relative population of the 
ground and first excited motional state in 1 5'o and thus 
the atomic ensemble temperature [l8l [l9j . As shown in 
Fig. 3(a), this measure yields an atomic temperature of 
5 /iK, in good agreement with time of flight tempera- 



ture measurements. Since the lattice trap depth is only 
~20 /iK, the harmonic approximation deteriorates for ex- 
cited motional states (n > 1, where n is the motional 
quantum number). One consequence of the anharmonic- 
ity of the optical confinement is the decreasing energy 
spacing between different motional states, spreading the 
inelastic sidebands into the individual n — > n + 1 tran- 
sitions. This anharmonic spreading partly explains the 
broad sidebands shown in Fig. 3(a). By knowing the trap 
depth and the energy of the n=0 — > n'=l transition, we 
solve for the exact eigenenergies of the cos 2 (/cz) depen- 
dent longitudinal lattice potential [2(J , where k is the lat- 
tice wave vector. The resulting eigenenergies, in units of 
lattice photon recoil energy, are shown in Fig. 3(b) as the 
lattice band structure. Knowing the relative frequency of 
each n — > n + 1 transition, we fit the blue sideband of Fig. 
3(a) with multiple Lorentzians to yield the amplitudes 
(assumed to be Boltzmann-distributed) and the widths 
of these transitions. The relative amplitudes among the 
various n — > n + 1 blue sidebands offer a more accurate 
measure of the atomic population of various longitudi- 
nal motional states and hence the sample temperature. 
This result, also 5 /iK, is consistent with measurements 
based on the relative blue and red sideband amplitudes. 
The individual n — > n + 1 transition linewidths are very 
broad compared to the elastic n — > n carrier transition. 
Although some intrinsic broadening is introduced by the 
band structure of the motional eigenstates (<1 Hz for 
lower states, ~kHz for upper states), major linewidth 
contributions arise from externally induced dissipations 
among motional states and radial variations of Q. Pos- 
sible external dissipation includes motional transitions 
induced by lattice phase and intensity noise, two-body 
atomic collisions, and the blackbody radiation. Unlike 
typical ion trapping experiments where individual mo- 
tional states are indistinguishable due to deep harmonic 
trapping, the scenario described here permits detailed 
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FIG. 3: (color online) (a) Spectroscopic trace of saturated, 
elastic 1 So - 3 Po electronic transition and the inelastic mo- 
tional sidebands at 80 kHz detuning from the center transi- 
tion, (b) Lattice band structure derived from exact solution 
of cos 2 (kz) dependent longitudinal lattice potential, (c) Fit 
to the blue sideband of (a) where the sideband is composed 
of multiple n — > n+ 1 transitions. Here the trap depth is ^20 
fiK while the atomic sample temperature is 5 fiK. The blue 
vertical bars denote the position of each n — > n + 1 transition 
from (b). (d) Narrow 1 S'o - 3 Po elastic carrier spectrum. 



studies of level-specific dissipation dynamics. 

Spectroscopy of the 1 6'o - 3 Po elastic carrier transition 
of 87 Sr in an optical lattice has the potential to yield one 
of the highest optical resolutions ever measured [2l|. Of 
course, a number of potential relaxation mechanisms can 
limit the achievable coherence time. A narrow spectrum 
of the carrier transition is shown in Fig. 3(d). The < 
200 Hz full width at half maximum linewidth allows for 
relatively rapid averaging of multiple spectra to achieve 
statistical uncertainties of the transition center frequency 
below the 5 Hz level, typically limited by the hydrogen 
maser stability used for frequency counting. Our efforts 
to reduce the broadening mechanisms were made in sev- 
eral different areas. To eliminate linewidth broadening of 
the fiber-optic-transferred probe laser (~100 Hz) caused 
by fiber phase noise, we implemented fiber noise can- 
cellation to sub- Hz precision |22|. Lack of perfect over- 
lap between the probe and lattice lasers can introduce 
Doppler broadening along the weakly confined transverse 
axes of the 1-D lattice. For this reason, we ensured that 
the well overlapped lattice and probe lasers had strong 
back coupling into both fibers from the retroreflecting 
lattice mirror. Additionally, the finite size of the probe 
beam creates a Heisenberg limited spread of k p which 



TABLE I: Typical systematic corrections and their associated 
uncertainties for the absolute frequency of the 1 So - 3 Po clock 
transition. 



Contributor 


Correction (Hz) 


Uncertainty (Hz) 


AC Stark shift (lattice) 


-17 


8.3 


AC Stark shift (probe) 


7.0 


0.9 


Blackbody shift 


2.2 


0.02 


Zeeman shift 





12 


Recoil/Doppler shift 





1.9 


Density shift 


-2.0 


13 


Gravitational shift 


-0.66 


0.14 


Cs-maser calibration 


46.8 


0.5 


Systematic total 


36 


20 



samples some of the transverse motion. This effect is be- 
low 50 Hz. In future work, these broadening mechanisms 
will be circumvented by implementation of a 3-D opti- 
cal lattice. State preparation of the atomic population 
in ground state mp sublevels was not attempted, mak- 
ing inhomogenous Zeeman broadening possible because 
of the differential shift of ~100 Hz/G x mp between 1 Sq 
and 3 Po states due to the latter's hyperfine mixing [23^ . 
We used Helmholtz coil pairs to minimize residual mag- 
netic fields. If not operated at A mag i C , the lattice laser can 
introduce large inhomogenous AC Stark broadening. By 
varying the lattice power and wavelength around A mag i C , 
we ensured linewidth contributions were below 30 Hz. 
Also, dissipation processes in the higher motional states 
can cause broadening in the elastic n — > n transition. By 
sideband cooling the population to the ground motional 
state, this effect could be reduced 24]. To date, we have 
not observed any effect of atomic density in the lattice 
on the measured linewidth. 

Figure 4 summarizes the measurement of the 1 So - 3 Po 
transition frequency within three months. Limited by a 
statistical uncertainty of 2.8 Hz, we have performed a de- 
tailed study of systematic shifts and their corresponding 
uncertainties, which are listed in Table 1. The most basic 
systematic error associated with spectroscopy in an op- 
tical lattice is the AC Stark shift of the clock transition 
due to the confining optical potential. To experimentally 
determine the magnitude and uncertainty of this shift, 
we varied the lattice intensity around our typical operat- 
ing intensity of Iq = 35 kW/cm 2 at wavelengths below, 
above, and near A mag i C . At our typical operating wave- 
length of 813.437 nm, Fig. 4(a) shows the frequency shift 
as a function of intensity. The slope of this shift yields 
an overall correction of 17(8.3) Hz at the typical lattice 
intensity of Jo- We have also experimentally determined 
that for the lattice intensity of io> the induced frequency 
shift is ~2 mHz for a lattice frequency deviation of 1 
MHz from the magic wavelength. Combining this with 
our measurement of the lattice AC stark shift of 17 Hz at 
813.437 nm yields A mag i C = 813.418(10) nm, in agreement 
with 0. 

With potentially high atomic densities in an optical 
lattice, characterization of density shifts is important. 
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FIG. 4: (color online) The measured 1 So 
quency versus (a) lattice intensity (Jo : 
atomic density, and (c) magnetic field, (d) JILA measure- 
ments over a 3 month period, with each data point represent- 
ing an averaged daily frequency measurement. The results 
reported in this work (lower bars) and in Ref. 6] (upper 
bars) are both shown with the total (outer box) and statisti- 
cal (inner shaded area) errors. 



We have ~400 lattice sites with typically ~250 atoms 
per site, yielding densities of ~10 12 atoms/cm 3 . Vary- 
ing the atomic density by a factor of 50 (Fig. 4(b)), we 
find the density shift at our typical operating density of 
2 x 10 12 atoms/cm 3 to be 2(13) Hz. Asymmetric popula- 
tion distributions among 1 S U ground state mp sublevels 
can lead to Zeeman shifts of the transition frequency. We 
measured a shift of 32 Hz/G (Fig. 4(c)). By keeping the 
magnetic field <400 mG during spectroscopy, the Zee- 
man shift uncertainty is 12 Hz. The probe beam itself 
can induce a frequency shift through two different phys- 
ical mechanisms. The first is the AC Stark shift of the 
1 Sq and 3 Pq levels due to their couplings to other states 



by the probe laser. For our investigation, the AC Stark 
shift due to the probe laser was exacerbated by using an 
electro-optic (EO) modulator to probe the atoms with 
a weak sideband while retaining an off resonant carrier 
for use in our fiber noise cancellation signal. The second 
frequency shift mechanism arises from the posssibility of 
a small probe beam misalignment with respect to the 
lattice laser, permitting photon recoil shifts in the trans- 
verse direction of the optical trap. We could separate 
these power-dependent effects by either varying the total 
698 nm light intensity incident on the atoms or changing 
the relative carrier-sideband amplitude through the EO 
modulation index. We determined each of these effects 
within 2 Hz uncertainty. 

The frequency reference used for the optical measure- 
ment was a hydrogen maser directly calibrated by the 
NIST Fl Cs fountain clock, available to us by an opti- 
cal fiber link from NIST to JILA |2£[. The approximate 
14 m lower elevation of our JILA Sr experiment relative 
to the NIST fountain clock introduced a < 1 Hz gravi- 
tational shift. Including all systematic effects discussed 
above, the total uncertainty (added in quadrature) is 20 
Hz. The final number we report for the 87 Sr 1 So - 3 Po 
transition frequency is 429, 228, 004, 229, 867 ± 20(sys) 
± 2.8(stat) Hz. We note this result disagrees by three 
standard deviations with that of Ref. 6] measured with 
a GPS-calibrated reference. 

We have presented the most precise spectroscopic mea- 
surement of the 87 Sr 1 S'o - 3 Po transition frequency in an 
optical lattice with extensive studies of systematic un- 
certainties. We apply this high resolution spectroscopy 
to investigating motional properties of lattice trapped Sr 
atoms. This work demonstrates the strength of the sys- 
tem for an optical atomic clock. Future work will improve 
further upon the reported precision and accuracy. 
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